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Abstract 



A fast and simple O(logn) iteration algorithm for individual Lucas numbers 
is given. This is faster than using Fibonacci based methods because of the 
structure of Lucas numbers. Using a \/5 conversion factor gives a faster Fi- 
bonacci algorithm because the speed up proposed in [5] also directly applies. 

A fast simple recursive algorithm for individual Lucas numbers is given 
that is O(logn). 
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t-h ; 

When initial conditions for the Fibonacci recurrence, F n = F n _i + F n _ 2 , 
are changed to 1,3 the result is the Lucas numbers (1891), L n = L n _i + L n _ 2 , 
with Li = 1, L 2 = 3. It is convenient to define L Q = 2, then L 2 = Li + L . 
There are a number of iterative solutions for this recursive definition, us- 
c5 1 ing addition, that are 0(n) [6], the lower bound for enumerating the sequence; 

however, individual numbers can be found in sublinear time. DeMoivre pub- 
lished a closed formula in 1730 that requires n multiplications [1]. F n was 
first shown to be found in 0(logn) time in 1978 [3 ] followed by improved 
algorithms [4,5]. None of these dealt directly with Lucas numbers perhaps 
because L n = F n _i + F n+ i, requiring two calls. 

For simplicity, we need not discuss the size of L n , which grows as 0(F n ), 
because this is the same for all algorithms (when the bit model executions 
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costs are excluded [4,5]). Execution costs are based on the represention 
operations, such as number of additions and multiplications. 

The basic idea for obtaining an O(logn) algorithm is the use of doubling 
formulas. Simple doubling only gives values of n that are powers of 2. Some 
manipulation is required to to solve for all n. Using the well known iden- 
tities [1], F n+m = F m * F n+1 + F m _ x * F n and F n+1 * F n ^ = F n 2 + (-1)" , 
algorithms can be developed for F n . A clever approach is to combine the 
related Fibonacci and Lucas sequences as for example, Fm = F n * L n and 
L 2n = L\ — 2, where 2n is even [3,4,5]. A single formula can only be used to 
find L n when n is a power of 2 and some manipulation is required to obtain 
a general solution for all n. 

Noting that the general doubling formula for Lucas is actually L 2 k = L\ — 
(— 1) * 2 and substituting k+1 for k then gives L 2 k+ 2 = L 2 k+1 + (—l) *2. This 
single fundamental equation form is then the basis for a compact solution 
compared to [5]. 

2. Iteration 

To calculate L n for any n requires three adjacent terms in the sequence. 
Using a method based on a recursive based calculation the two terms to 
double are selected by markOdd() in the following algorithm. 

2.1. Middle 

For simplicity, the sequence variables are renamed as follows: LL = 
L 2 k,LM = L 2 k+i,LH = L 2 k+ 2 The following formulas are used to double 
until L n is found. 

LL= LL*LL - p*2 
LH= LM*LM + p*2 
LM= LH-LL 

Depending on the step, two different updates are required to prepare for 
the next step these are selected using an odd/even tag. This can be seen in 
the algorithm. 

Two square multiplications, two sign calculations, and three additions in 
the iterative loop make this very fast compared to [5]. No pre conditions nor 
post conditions are used as was the case in [5]. All general methods need 
some kind of odd(i) selected calculation, when n is not a power of 2, and 
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this has been reduced here to the insertion of a shift forward of two sequence 
terms and one assignment. 

The algorithm can be expressed as a recursion in a direct way, which 
is done in the next section. Using the relation, F n = \(L n /\fE — 0.5) , our 
algorithm also gives a very fast and simple calculation of a Fibonacci number. 

Algorithm : Middle(n) given n > 1, return LL = L n 
N 4- [lgn 

array markOdd[N] <— 
i <— n; j <— N 

while(j > 0){if(odd(i)) markOdd(j) — l;i — i/2;j = j — 1} 

whileQ <N-1) 
p 4- 1 

if(markOdd(j)){LL <- LM;LM <- LH;p <- -1 } 
LL ^— LL * LL — p * 2 
/.// <- LM * LM + p*2 
LM <- LH - LL 

3^3 + i; 
endwhile 

if(markOdd(j)) LL <- LM 
return LL 

3. Recursion 

The recursion L 2n = 1^ — 2 when n is even can be used to find L n 
in O(logn) time when n is a power of 2. This can be executed recursively 
by using L k = L 2 k ^ 2 — 2 The general form of this doubling formula for 

Lucas is L 2 k = L\ — (— l) k * 2 [1] and substituting k+1 for k then gives 
L 2 k+2 = L\+\ + {— l) fc * 2. The result on two adjacent position terms is two 
even position semi-adjacent terms. To obtain L n when n is odd, calculate 
L n+ i — L n _i. Thus to calculate L n for any n requires three adjacent terms 
in the sequence. 

3.1. Ripple 

For simplicity, the sequence variables are renamed as follows: LL = 
Lk, LH = Lfc+2 The following formulas are used to double until L n is found. 
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LL= LL*LL - p*2 
LH= LM*LM + p*2 



Depending on the odd and even values of n and n/2, different updates 
are required. The correct equation to use is selected by the variable p, as can 
be seen in the algorithm. 

Algorithm : Ripple{n) given n > 1, 

return L n 

if(n = 2) return 3 

if(n = 3) return 4 

if(n = 4) return 7 

if (even(n/2)) then p 1 else pi f 

LL <- Ripple([n/2\) 

if (even(n)) then return LL * LL — p * 2 

else LH <- Ripple(\n/2]) 

return LH * LH — LL * LL + p * 4 

If a compiler is able to remember multiple identical calls, the algoritm 
can be simplified to remove LH and HL. 

if (even(n)) then return L{n/2) 2 — p * 2 
else return L( [n/2] ) 2 - L( \n/2\) 2 + p * 4 

Each call uses a selection of p, one or two square multiplications, and 
one signed addition. Making this very fast compared to [5] , when the cost of 
recursion is ignored. No pre conditions nor post conditions are used as was 
the case in [5], except to define the recursion basis. 

4. Conclusions 

For linear Fibonacci algorithms, a change in initial conditions gives a 
Lucas algorithm. Any O(logn) Fibonacci algorithm can be used to find 
the Lucas number indirectly by the relation L n = F n _i + F n+1 . While 
this requires two external calls, these algorithms can usually be modified 
internally to produce L n in a single call (which may only hide the two external 
calls). Conversely, Fibonacci numbers can be derived from Lucas numbers 
using a divison factor. By the nature of the simple foundation relation 
of Lucas squares, our algorithm is simpler and faster than other Fibonacci 
algorithms. 
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Since our algorithm only uses two square multiplication per iteration, the 
speed argument for squares in [5], claimed to be the fastest for very large n, 
also applies to our method. This algorithm is much simpler and faster by a 
constant factor than that in [5] (which did not internally implement the fast 
FFT-based Schoenhage-Strassen process for fast multiplication). Thus, the 
two algorithms can be directly compared. 

In Ripple the the average number of multiplications is less than 2. The 
algorithm only uses squares, and the argument in [5], claimed to be the fastest 
for very large n, also applies to this method. 

While Lucas numbers can be derived from finding Fibonacci numbers, it 
is always more efficient to calculate the Lucas numbers directly. 
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